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ABSTRACT 

We use three samples (3CRR, 6CE and 6C*) selected at low radio frequency to con- 
strain the cosmic evolution in the radio luminosity function (RLF) for the 'most lu- 
minous' steep-spectrum radio sources. Although intrinsically rare, such sources give 
the largest possible baseline in redshift for the complete flux-density-limited samples 
currently available. Using parametric models to describe the RLF which incorporate 
distributions in radio spectral shape and linear size as well as the usual luminosity and 
redshift, we find that the data are consistent with a constant comoving space density 
between z ~ 2.5 and z ~ 4.5. We find this model is favoured over a model with similar 
evolutionary behaviour to that of optically-selected quasars (i.e. a roughly Gaussian 
distribution in redshift) with a probability ratio of ~ 25 : 1 and ~ 100 : 1 for spatially- 
flat cosmologies with = and f^A = 0.7 respectively. Within the uncertainties, 
this evolutionary behaviour may be reconciled with the shallow decline preferred for 
the comoving space density of flat-spectrum sources by Dunlop & Peacock (1990) and 
Jarvis & Rawlings (2000), in line with the expectations of Unified Schemes. 

Key words: galaxies:active - galaxies: luminosity function, mass function - radio 
continuum : galaxies 



1 INTRODUCTION 

Speculation concerning the cosmic evolution of the ra- 
dio source population has been at the forefront of as- 
trophysical research since the 1960s. It soon became 
clear that the radio luminosity function (RLF) declines 
steeply towards high radio luminosities, and that the 
(comoving) space densities of the rarest, most luminous 
quasars and radio galaxies were much higher at epochs 
corresponding to z ~ 2 than they are now (Longair 
1966). As emphasised by Peacock (1985), Dunlop & 
Peacock (1990; hereafter DP90), Shaver et al. (1996; 
hereafter SH96), Jarvis & Rawlings (2000; hereafter 
JROO) and many others the crucial advantage of any 
radio-based work is that with sufficient optical follow- 
up, it can be made free of optical selection effects, such 
as increasing dust obscuration at high redshift. 

The term 'redshift cut-off made an early appear- 
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ance in the literature (e.g. Sandage 1972), and is now 
taken to mean any significant decline in the space den- 
sity over 2.5<z<5. Peacock (1985) and then DP90, 
found evidence for a redshift cut-off in the distribution 
of flat-spectrum radio quasars over the redshift range 
2 — 4. Through failing to find any flat-spectrum radio 
quasars at z > 5 in a survey covering ~4sr, SH96 (see 
also Shaver et al. 1998; hereafter SH98) argued for a 
drop in space density between 2.5 6 of more than 
1 dex. JROO re-addressed this question using a similar 
sample to the flat-spectrum sample of SH96/SH98. By 
considering the effects of a distribution in spectral in- 
dex and also the presence of any curvature in the ra- 
dio spectra, JROO found evidence for a decline between 
2.5 < 2; < 5 by a factor 4 in agreement with DP90, 
although they found that steeper declines such as that 
preferred by SH96, as well as models invoking a con- 
stant comoving space density at high redshift, can only 
be ruled out at the 90% confidence level with the sur- 
veys currently available. 
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DP90 also claimed the first evidence for a simi- 
lar redshift cut-off in the steep-spectrum radio popula- 
tion. According to simple Unified Schemes (e.g. Barthel 
1989; Antonucci 1993) steep-spectrum radio sources are 
the parent, unbeamed, population for Doppler-boosted 
flat-spectrum objects, and should show similar evolu- 
tionary behaviour (Jackson & Wall, 1999)|. However, 
the high selection frequency (2.7 GHz) of the steep- 
spectrum samples used by DP90 meant that their sur- 
vey was sensitive to only the most luminous sources at 
high redshift, an effect which is made more pronounced 
by the positive correlations between spectral index and 
redshift /radio luminosity (e.g. Blundell, Rawlings & 
Willott 1999). The necessity of considering the effects 
of the distribution in radio spectral parameters, high- 
lighted by JROO for the flat-spectrum population, is es- 
pecially acute for this population which typically have 
steep and curved spectra. There are also worries about 
the incompleteness of the faint radio-selected sample 
used by DP90 (see their Appendix A). Moreover, DP90 
were forced to use photometric redshift estimates for 
most of their high-redshift steep-spectrum objects, and 
since these were based on the still poorly-understood 
K — z diagram for radio galaxies (e.g. Bales et al. 1997; 
Jarvis et al. 2001b), this leads to considerable uncer- 
tainty in their results. 

To ameliorate these problems, our group has been 
seeking spectroscopic redshifts for samples selected at 
low radio frequency from the 6C and 7C catalogues (see 
e.g. Rawlings et al. 1998). There are now three low- 
frequency-selected samples with essentially complete 
redshift information: 3CRR (Laing, Riley & Longair, 
1983); 6CE (Bales et al. 1985; Rawlings, Bales & Lacy, 
2001) and the 7C redshift survey (7CRS; BlundeU et 
al. in prep.; Willott et al. in prep.). Investigation of the 
RLE of steep-spectrum quasars (Willott et al. 1998), the 
first of its kind, found no definitive evidence for a red- 
shift cut-off, albeit with only a small sample of quasars. 
Using the combined 3CRR, 6CB and 7CRS dataset, 
comprising 357 sources, Willott et al. (2001) obtained 
similar results for the entire low-frequency population. 

In this paper we focus on whether there is any ev- 
idence for a redshift cut-off in the space density of the 
most radio-luminous steep-spectrum radio sources. Our 
work on follow-up spectroscopy of the 6CB (Rawlings 
et al. 2001) and 6C* (BlundeU et al. 1998; Jarvis et 
al. 2001a, 2001b) samples provides redshifts for well 
defined samples containing the distant counterparts of 
the most luminous radio sources in the 3CRR survey 
(Laing, Riley & Longair 1983). This analysis of the RLP 
is different from that by Willott et al. (2001) since it 
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(GPS) sources also play a major role in the flat-spectrum surveys 
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also includes 6C*, which is crucial in sampling to high 
redshift. 

In Sec. H we summarise the relevant properties of 
the 3CRR, 6CB and 6C* samples, and describe how 
the most luminous radio sources have been drawn from 
them. In Sec. ^ we highlight the reasons for concen- 
trating on just the most radio-luminous sources and 
in Sec. ^ we describe the parametric models used in 
a maximum likelihood analysis which yield estimates 
of the RLF. We use the banded V/Vmax statistic in 
Sec. H to compare our parametric modelling with a non- 
parametric approach and in Sec. ^ we compare our pre- 
dicted redshift distributions with the data and highlight 
the importance of small number statistics. In Sec. we 
discuss our results in the context of unified schemes of 
radio AGN. 

We assume Ho — 50 km s^^ Mpc^^, and obtain 
all results in two spatially-flat cosmologies (J^m = 1 
and f2A — 0, hereafter cosmology I and f2M = 0.3 and 
f^A = 0.7, hereafter cosmology II). We use the conven- 
tion for spectral index a that Si, cx y~°' and aJ^J rep- 
resents the spectral index between observed frequencies 
vi and V2 ■ All luminosities quoted are measured in units 
of WHz^^sr~^at rest-frame 151 MHz unless stated oth- 
erwise. 



2 COMPLETE SAMPLES OF RADIO 
SOURCES SELECTED AT LOW 
FREQUENCY 

2.1 The advantages of using three samples 

Flux-density-limited samples over a particular sky area 
provide a direct way of constraining the RLF if they 
have complete redshift information. In this analysis, in 
which only the 'most luminous' sources are considered, 
the 7CRS with a flux-density limit 5'i5i > 0.5 Jy cov- 
ering a total sky area of 0.022 sr would contribute only 
one such source, principally due to the small sky area 
it surveyed and the steepness of the RLF (see Willott 
et al. 2001). Thus, inclusion of 7CRS was deemed un- 
necessary for the purpose of this paper and we only use 
three samples, 3CRR, 6CB and 6C* which we describe 
in the next three sections. 

2.2 3CRR 

We use the 3CRR sample of Laing, Riley & Longair 
(1983), which contains 173 sources. Following Blundell 
et al. (1999) we exclude 3C231 (M82) since the radio 
emission is a consequence of star-formation (Condon 
et al. 1990) rather than an active nucleus. Two other 
sources are also excluded (3C 345 and 3C 454.3) because 
their emission only exceeds the 178 MHz flux-density 
limit because of Doppler boosting of their core emission. 
The sample has 100% spectroscopic redshifts. 
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Sample 


Flux-density 


Sky Area 


Number of 


Uncertain 


Spectroscopic 


Reference 




limits / Jy 


/ sr 


Sources 


redshifts 


completeness 




3CRR 


5l78 > 10.9 


4.24 


170 





100% 


Laing et al. (1983) 


6CE 


2.00 < 5151 < 3.93 


0.10 


58 


7 


97% 


Rawlings et al. (2001) 


6C* 


0.96 < 5i5i < 2.00 


0.13 


29 


6 


100% 


Jarvis et al. (2001a) 



Table 1: Summary of the three samples used in the analysis to constrain the RLF of the most luminous radio 
sources. The spectroscopic completeness figure neglects uncertainties in a small number of redshifts in the 6CE and 
6C* samples. 
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Rodshift. 

Figure 1: Rest-frame 151 MHz luminosity (Lisi) ver- 
sus redshift z plane for the three samples used in our 
analysis in cosmology I; 3CRR (circles), 6CE (triangles) 
and 6C* (stars). The rest-frame 151 MHz luminosity 
Li5i has been calculated according to the polynomial 
fit to the radio spectrum described in Sec. 4.2. The 
6C* source 6C*0140-h326 at z = 4.41 is represented 
by two points joined by a solid vertical line. This is 
because 60*0140-1-326 is thought to be subject to grav- 
itational lensing by a foreground galaxy (Rawlings et 
al. 1996) and the highest luminosity point corresponds 
to the lensed luminosity whereas the lower point corre- 
sponds to the de-magnified, intrinsic luminosity of the 
source in the absence of significant lensing. The area be- 
tween the horizontal lines is the region which contains 
the 'most luminous' sources according to our definition. 
The curved lines show the lower flux-density limit for 
the 3CRR (dotted) and 7CRS (solid). The dashed lines 
correspond to the limits for the 6CE sample and the 
shaded region shows the 6C* flux-density limits (all as- 
suming a radio spectral index of 0.5). Note the area be- 
tween the 3CRR sources and 6CE sources contains no 
sources, this is the area which corresponds to the ab- 
sence of a flux-density-limited sample between the 6CE 
(S'lsi < 3.93 Jy) and 3CRR (5i78 > 10.9 Jy) samples. 



2.3 6CE 

The 6CE sample originally defined by Eales et al. 
(1985), and recently revised by Rawlings et al. (2001), 
probes flux-densities a factor ^ 6 fainter than 3CRR. 
The sample consists of 59 objects, however, one of these 
(6CE1036-I-3616) is excluded for the statistically unbi- 
ased reason of there being a bright star occluding any 
likely optical identiflcation. The sample considered here 
consists of 58 objects, with only two objects lacking a 
spectroscopic redshift, a completeness of ^ 97%, with a 
further 7 redshifts not yet secure. 

2.4 6C* 

The 6C* sample (Blundell et al. 1998) is a filtered sam- 
ple from part VI of the 6C survey of radio sources at 
151 MHz which was designed to find radio sources at 
z > 4. The completeness of the sample is summarised 
in Table 0. At the faint flux-density limit of 6C*, the 
source counts of radio sources at 151 MHz are of or- 
der 10^ per steradian, making it very difficult to con- 
struct a complete sample with spectroscopically deter- 
mined redshifts over large solid angles. Moreover, the 
low-luminosity and thus low-redshift sources supply a 
significant contribution to the source counts at these 
flux-density limits making it increasingly difficult to use 
faint flux-density-limited samples to probe evolution of 
radio sources at high redshift. To combat this prob- 
lem, additional selection criteria based on radio prop- 
erties can be used to filter out the low-luminosity (low- 
redshift) sources. These additional selection criteria are 
discussed in detail in Bhmdell et al. (1998), but to sum- 
marise, only sources with spectral indices steeper than 
o^aisK? — 0.981, and angular size 6 < 15arcsec are 
included. These selection criteria introduce additional 
complications in any analysis where the 6C* sample is 
used, and need to be treated with careful considera- 
tion as will be discussed in Sec. |4.1| . This sample has 
now been reduced from 34 to 29 objects since images 
from the NVSS survey (Condon et al. 1998) have con- 
firmed the suggestions of Blundell et al. (1998) that 
these objects were the hotspots of large {9 > 15") dou- 
ble sources. The optical spectroscopic follow-up of this 
sample has now been completed (Jarvis et al. 2001a) 
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and of the remaining objects, 23 have secure spectro- 
scopic redshifts, and 6 have tentative redshifts. The red- 
shifts of these six sources are typicaUy based on a single 
definite emission hne identified with a specific feature 
on the basis of a rcdshift estimate from the if— band 
magnitude. 

It is also worth noting that the absence of strong 
Lya emission in the objects from 6CE and 6C* which 
lack secure redshifts implies that these source are prob- 
ably below z < 1.8, and should not affect the analysis 
in this paper (as we are only considering the most lu- 
minous sources which correspond to high redshifts at 
these flux-densities). 



in any determination of how the comoving space den- 
sity of radio sources evolves, we need to account for the 
sources excluded at all redshifts. This is crucial because 
if a significant number of sources between 2 < z < 3 
arc filtered out of the sample, and objects at z > 3 
are not, then on first inspection the results would be 
very misleading as the derived comoving space density 
at high redshift would be higher relative to the comov- 
ing space density at z ~ 2. To incorporate the effects 
of the selection criteria into the parametric modelling 
we adopt functional forms which parameterise the dis- 
tributions in spectral shape and linear size, in addition 
to luminosity and redshift. These functional forms are 
described in the following sections. 



3 SELECTING 'THE MOST LUMINOUS 
RADIO SOURCES' 

The principal aim of this paper is to investigate the evo- 
lution in the comoving space density of steep-spectrum 
radio sources directly. Thus, we concentrate on only the 
'most luminous' sources, which we are able to detect to 
z > 4 in our flux-density-limited samples. Another rea- 
son for concentrating on just the most luminous objects 
is that in such a narrow luminosity range any bias from 
intrinsic correlations which are present in radio samples, 
such as luminosity - spectral index and linear size - spec- 
tral index correlations (e.g. Blundell et al. 1999) will be 
minimised. Selecting which sources from the three sam- 
ples should be included in our analysis as the 'most 
luminous' was simply a question of where to place the 
lower limit in radio luminosity. To enable us to com- 
pare results between various cosmologies we are forced 
to alter the limits between cosmologies. Therefore to 
make the analyses consistent we only consider the top 
decade in L151 in cosmology I [VLm = 1 and 17a ~ 0). 
Thus all sources with logj^Q L151 > 27.63 were included. 
For cosmology II {Q,m — 0.3 and SIa = 0.7) we redefine 
our lower limit in luminosity at logj^Q Lisi > 27.90, this 
is not strictly the lower limit of the top-decade (which 
is logio ^151 — 27.95) but allows the same number of 
sources into the modelling process as in cosmology I. 

This selection process led to the inclusion of 57 
3CRR sources for cosmology I and 55 3CRR sources for 
cosmology II, four (five) 6CE sources, and five (six) 6C* 
sources in cosmology I (cosmology II). These sources are 
tabulated in Table ^. 



4 MODELLING THE RLE 

4.1 6C* - what are we missing? 

One of the principal problems in dealing with a filtered 
sample such as 6C* is determining which populations 
are filtered out. One of the aims of the 6C* sample was 
to find objects at z > 4. However, to use this sample 



4.2 Radio spectral shape 

Many radio spectra are known to exhibit curvature (e.g. 
Laing & Peacock 1980; Blundell et al. 1999) thus ren- 
dering the use of a power-law spectral index inadequate 
in describing the radio spectra over a large frequency 
baseline. As a consequence of this, determining the rest- 
frame luminosity becomes less straightforward. Thus, to 
account for spectral curvature we fit the radio spectra 
with a polynomial of the form 

N 

y = logio 5"^ = X!"*^*' f^-^) 

where x = logio(j^ / MHz). We fitted the flux-density 
data using a Bayesian polynomial regression analysis 
which assessed the posterior probability density func- 
tion (pdf) for the required order of polynomial fit (e.g. 
Sivia 1996). This is the same method as used in Blun- 
dell et al. (1998) and Blundell et al. (1999). The results 
of this analysis are presented in Table H; note that in 
no cases was N > 2 preferred, so the fits are either 1st 
or 2nd order polynomials. 

4.3 Parametric modelling of the RLE for the 
most luminous radio galaxies and quasars 

To extract as much information as possible from our 
dataset, we use three parametric models, with as few 
free parameters as possible, to estimate the RLF of the 
most luminous low-frequency-selected radio sources. 

Following JROO we construct a luminosity function 
separable in luminosity and redshift along with addi- 
tional parametcrisations describing the distributions in 
spectral shape and projected linear size. 

We use a single power-law to parameterise the lu- 
minosity function pl{Lizi), i.e. 

PL(ii5i)=(^) (2) 

where /3 is a dimensionless free parameter and Lo is 
a normalising luminosity fixed at the lower limit of the 



The steep-spectrum redshift cut-off 5 



loglo(-^151 / W Hz sr ) ai a2 D j kpc Optical 

Cos I Cos II Cos I Cos II type 



3C9 


2 


012 


28 


63 


28 


94 


-0.47 


-0 


10 


114 


5 


164 





Q 


3C13 


1 


351 


28 


06 


28 


33 


-0.33 


-0 


11 


241 


7 


330 


8 


G 


3C14 


1 


469 


28 


07 


28 


35 


-0.44 


-0 


08 


205 


5 


284 


1 


Q 


3C22 





938 


27 


73 


27 


96 


-0.58 


-0 


05 


206 


2 


269 


2 


B 


3C43 


1 


457 


28 


09 


28 


37 


-0.75 





00 


21 


4 


29 


6 


Q 


3C55 





734 


27 


75 


27 


95 


-0.09 


-0 


16 


572 


9 


723 


8 


G 


3C65 


1 


176 


27 


94 


28 


19 


0.16 


-0 


18 


149 


7 


201 


4 


S 


3C68.1 


1 


238 


28 


00 


28 


26 


-0.36 


-0 


10 


447 


9 


606 


4 


Q 


3C68.2 


1 


575 


28 


23 


28 


52 


-0.39 


-0 


16 


189 


7 


264 


5 


G 


3C147 





545 


27 


77 


27 


94 


0.76 


-0 


22 


22 


1 


26 


8 


Q 


3C175 





768 


27 


77 


27 


97 


-0.92 


-0 


01 


391 


3 


497 


4 


Q 


3C181 


1 


382 


28 


12 


28 


39 


-0.23 


-0 


11 


49 





67 


2 


Q 


3C184 





994 


27 


76 


28 


00 


0.03 


-0 


16 


40 


8 


53 


7 


G 


3C186 


1 


063 


27 


89 


28 


14 


0.18 


-0 


22 


13 


7 


18 


2 


Q 


3C190 


1 


197 


28 


04 


28 


30 


-0.63 


-0 


04 


57 


7 


77 


8 


Q 


3C191 


1 


952 


28 


50 


28 


81 


-0.78 


-0 


03 


40 


3 


57 


6 


Q 


3C196 





871 


28 


35 


28 


58 


-0.07 


-0 


13 


83 


6 


108 





Q 


3C204 


1 


112 


27 


94 


28 


19 


-1.18 





01 


314 


2 


419 


6 


Q 


3C205 


1 


534 


28 


19 


28 


48 


-0.34 


-0 


11 


153 


5 


213 


4 


Q 


3C208 


1 


109 


28 


03 


28 


28 


-0.26 


-0 


13 


94 


4 


126 


1 


Q 


3C212 


1 


049 


27 


92 


28 


16 


-0.78 





00 


77 





102 





Q 


3C217 





897 


27 


66 




— 


-0.33 


-0 


11 


100 


8 






G 


3C226 





817 


27 


72 


27 


94 


-0.52 


-0 


08 


289 


1 


370 


6 


G 


3C239 


1 


781 


28 


41 


28 


71 


-0.48 


-0 


10 


93 


7 


132 


5 


G 


3C241 


1 


617 


28 


11 


28 


40 


0.36 


-0 


24 


7 


7 


10 


8 


G 


3C245 


1 


029 


27 


84 


28 


08 


-0.67 





00 


77 


7 


102 


7 


Q 


3C252 


1 


105 


27 


95 


28 


20 


-1.03 


-0 


01 


514 


9 


687 


1 


G 


3C254 





734 


27 


72 


27 


92 


-0.20 


-0 


14 


105 


7 


133 


5 


Q 


3C263.1 





824 


27 


76 


27 


98 


-0.21 


-0 


12 


56 


2 


72 


2 


G 


3C265 





810 


27 


86 


28 


07 


-0.82 


-0 


03 


643 


3 


823 


6 


S 


3C266 


1 


272 


27 


97 


28 


24 


-0.34 


-0 


11 


38 


8 


52 


7 


G 


3C267 


1 


144 


28 


00 


28 


25 


-0.60 


-0 


05 


326 


7 


437 


9 


G 


3C268.1 





973 


27 


97 


28 


21 


-0.70 





00 


390 


5 


512 


4 


G 


3C268.4 


1 


400 


27 


98 


28 


26 


-0.30 


-0 


10 


93 


6 


128 


7 


Q 


3C270.1 


1 


519 


28 


24 


28 


53 


-0.96 





02 


102 


5 


142 


2 


Q 


3C280 





996 


28 


05 


28 


29 


-0.50 


-0 


05 


123 


4 


162 


4 


G 


3C287 


1 


055 


27 


76 


28 


00 


0.27 


-0 


13 





8 


1 





Q 


3C289 





967 


27 


71 


27 


94 


-0.14 


-0 


12 


84 


8 


111 


2 


G 


3C294 


1 


786 


28 


38 


28 


68 


-0.82 


-0 


05 


125 


5 


177 


4 


G 


3C295 





461 


27 


81 


27 


96 


0.79 


-0 


26 


37 


7 


44 


9 


G 


3C309.1 





904 


27 


83 


28 


06 


-0.00 


-0 


10 


24 


4 


31 


7 


Q 


3C318 


1 


574 


28 


11 


28 


40 


-0.01 


-0 


14 


6 


8 


9 


5 


Q 


3C322 


1 


681 


28 


20 


28 


49 


-0.54 


-0 


07 


278 


6 


391 


2 


G 


3C324 


1 


206 


28 


05 


28 


31 


-0.15 


-0 


14 


86 


1 


116 


2 


S 


3C325 





860 


27 


72 


27 


95 


-0.26 


-0 


10 


133 


4 


172 


2 


G 


3C336 





927 


27 


71 


27 


94 


-0.68 


-0 


03 


183 


1 


238 


7 


Q 


3C356 




























Q 


4C16.49 


1 


296 


28 


05 


28 


31 


-1.00 





00 


137 


8 


187 


6 


Q 


4C13.66 


1 


450 


28 


10 


28 


38 


-0.81 





00 


51 


4 


71 





G 


3C368 


1 


132 


28 


06 


28 


31 


-0.50 


-0 


13 


67 


9 


90 


9 


G 


3C380 





691 


28 


09 


28 


29 


-0.45 


-0 


04 


158 


8 


199 





Q 


3C427.1 





572 


27 


64 




— 


-0.57 


-0 


07 


209 


6 




_ 


G 


3C432 


1 


805 


28 


31 


28 


61 


-0.38 


-0 


11 


108 


6 


153 


7 


Q 


3C437 


1 


480 


28 


12 


28 


40 


-0.01 


-0 


13 


294 


3 


407 


4 


G 


3C454 


1 


757 


28 


33 


28 


63 


-0.90 





00 


10 


9 


15 


4 


Q 


3C469.1 


1 


336 


28 


14 


28 


42 


-1.12 





00 


636 


7 


870 


3 


G 


3C470.0 


1 


653 


28 
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Table 2: Sources selected from the 3CRR, 6CE and 6C* samples which correspond to our definition of the 'most 
radio luminous'. The rest-frame 151 MHz luminosities were calculated by fitting a 1st or 2nd order polynomial to 
the available radio data as described in Sec. 4.2. The classes G, Q, B and S represent quasars, radio galaxies, broad- 
line radio galaxies and sources with scattered broad lines respectively, where we have followed the classifications of 
Willott et al. (1998) and Cos I and Cos II denote cosmologies I and II respectively. The solid horizontal lines are 
used where the source was not used for the specific cosmology due to it falling below our lower luminosity limit. 
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top-decade in luminosity, which is well above the 'break' 
luminosity in the RLF (Willott et al. 2001) and L151 is 
the rest-frame 151MHz luminosity. 
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Figure 2: The correlated parameters oi and 02 which 
describe the shape of the radio spectra. The circles 
are 3CRR sources; triangles correspond to 6CE sources 
and stars correspond to sources from the filtered 6C* 
sample. The inner (outer) contour corresponds to the 
la (2(t) spread about the correlation using the best-fit 
value of model B in cosmology I. 
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Figure 3: Histogram of the distribution in 
logig(Z)/kpc), where D is the projected linear 
size, for the top-decade in luminosity for cosmology I 
(top) and cosmology II (bottom). The white regions 
correspond to 3CRR sources, the light shaded regions 
are the contributions from the 6CE sample and the 
dark shaded regions are the 6C* sources. The solid line 
represents the Gaussian fit to the data from model A 



(Sec. 4.4) with an assumed normalisation factor which 



is incorporated within po in the parametric modelling. 

JROO highlighted the importance of a distribution 
in spectral index in their parametric modelling and also 



suggested that curvature would need to be incorporated 
in a full analysis of the RLF. The small number of 
sources in the JROO sample restricted the models to 
a small number of free parameters and the effects of 
spectral curvature could not be incorporated into their 
parametric models. However, we are able to introduce a 
distribution for spectral curvature into the modelling of 
the RLF here, due to the larger number of sources in the 
top-decade of luminosity from our set of low-frequency- 
selected samples. 

We fit the radio spectra with a polynomial of the 
form described in Sec. Fig. |2| shows that the pa- 
rameters describing spectral shape, oi and 02 are lin- 
early anti-correlated. Therefore to incorporate the ef- 
fects of spectral curvature into the modelling of the RLF 
we use a functional form which accounts for this anti- 
correlation and the spread about it. To parameterise oi 
and 02 we rotate the 01 — 02 plane by the free parameter 
ljj to transform the correlation so it becomes parallel to 
the fli axis, i.e. 

a'l - ^^a\ = [ai - f^a[)cos{uj) + (02 - fia'Jsm{uj), 

a'2 - Ma' = -(oi - fia' ) Sin(tj) + (02 - ) COS(W), (3) 



where ai and 02 are defined in Sec. 42, a'^ and a'2 are 
the transformed coordinates of ai and 02 rotated by the 
angle to and iia[ and /i^^ are the Gaussian peaks for the 
distributions in a'l and a'2 respectively and are invari- 
ant with the rotation. We can now fit the distributions 
in both a'l and a'2 with a Gaussian to characterise the 
spread about fitted mean values fia' and fia' : i-e. 



Pi (a'l) = cxp < - 




1 / «! - ^J■a[ 



(4) 



(5) 



where p,a'^ and fia[ are as defined previously, and CTq'^ 
and an' are the characteristic widths of the Gaussians 
describing the distributions in a'l and a'2 respectively. 
Thus, we obtain. 




Pa{ai,a2) = pi{a'i)p2{a'2). 



(6) 



To account for the angular size filtering in 6G* we 
also require a functional form which characterises the 
distribution in projected linear size D. For this we as- 
sume that there are no correlations between D and the 
other parameters in the small luminosity range we con- 
sider^. This enables us to model the linear size distri- 

t We note throughout this section that our chosen parameterisa- 
tions are imperfect, e.g. there are hkely to be weak correlations 
between radio spectral parameters, linear size and redshift even 
over a narrow range of radio luminosity (see BRW), but attempts 
to account for such subtleties would lead to an unacceptable in- 
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bution with a simple Gaussian distribution in log^^g 
i.e. 

where logj^Q D is the log of the projected linear size in 
kpc. logj^o Do is the peak of the Gaussian distribution 
and logj^o -^i is charaeteristie width of the Gaussian. 

The final factor needed to model the distribution of 
sources is a component describing the redshift distribu- 
tion of the most luminous radio sources. We use three 
independent models to achieve this. 

Model A is parameterised as a single Gaussian dis- 
tribution in redshift which is consistent with previous 
forms describing the evolution in the comoving space 
density from the literature (e.g. SH96), although as 
noted by JROO this form has the disadvantage of being 
a symmetrical function, thus coupling the low-redshift 
evolution with the high-redshift evolution for which 
there is not necessarily any physical basis. For model 
A we take 



Pa{z) = exp • 



Zl 



(8) 



where Zo is the redshift of the Gaussian peak and zi is 
the characteristic width of the Gaussian. 

Model B is a 1-tailed Gaussian which becomes con- 
stant beyond the Gaussian peak, this corresponds to 
a constant comoving space density above the Gaussian 
peak redshift, i.e. 



Pb{z) = 



exp 
1.0 



: z < Zo 
: z > Zo, 



(9) 



where Zo and zi are as defined previously. 

The third, and final model (C) uses a 1-tailed Gaus- 
sian to parameterise the low-redshift comoving space 
density and a power-law distribution at high redshift. 



I.e. 



Pc{z) = 



exp|-i(^ 



: z < Zo 
: z> Zo, 



(10) 



where Zo is the 'break' redshift where the model 
switches from the low- to the high-redshift form, zi 
is the characteristic width of the half-Gaussian and r] 
is the power-law exponent describing the high-redshift 
comoving space density. This model has the advantage 
of being free to carry on increasing at high redshift, as 
we assume there is no a priori reason for a decline over 
the redshift range of interest. However, it has the dis- 
advantage of an additional free parameter to the other 



models. Note that model B is just model C with ij fixed 
at zero. 

The distribution for the complete radio luminosity 
function p, is thus given by the product of all these 
factors, i.e. 

p{Li5i,z,ai,a2,D) = po x Pl(Li5i) x px{z) 

xpa{ai,a2) X pd{D), (11) 

where Po is the normalising factor and a free parame- 
ter measured in units of Mpc~^ and PL(-ti5i), Px{z) 
(for model X, where X G [A, B, C] ), Pa{ai, CL2) 
and Pd{D) are dimensionless distribution functions 
per (A log^Q iisi), per (A2), per (AaiAa2), and per 
(A log^o D) respectively. 

To find the best-fit values of the free parameters for 
the various models, the maximum likelihood method of 
Marshall et al. (1983) was used. If S is defined as — 21n£, 
where C is the likelihood function, then the aim is to 
minimise the value of S, which is given by: 

N 



S = -2y~^ln[p(Li5i., Zj, ai-, 02^, A 



mil 

xO(Li5i, z, oi, 02, D) 
dV 

X — d(logio Li5i ) dz dai da2 d(logio D) , 



(12) 



where p{Li5i, z, ai, a2, D) is the model distribution 
being tested, f7(Li5i, z, oi, 02, D) is the sky area avail- 
able for the samples under consideration in this paper 
for a given flux-density, and {AV/Az) is the differential 
comoving volume element. The first term is simply the 
sum over N som-ces in the defined sample. The sec- 
ond term is the integral of the model distribution being 
tested and should give w 2A'' for good fits. 

The lower and upper limits of the integral are 
27.63 < logio Li5i < 28.63 for cosmology I, and 27.90 < 
logio -^151 < 28.95 for cosmology II, < z < 10, -2.2 < 
ai < 1.0, -0.4 < 02 < 0.2 and -0.3 < log^o £> < 4.0. 
To illustrate the goodness-of-fit of the models the 2- 
dimensional Kolmogorov-Smirnov (KS) test, originally 
proposed by Peacock (1983) was used in the form of 
Press et al. (1992). The KS-test enables one to deter- 
mine the probability, PkSj that a model distribution 
is a true representation of the dataset. However, the 
2D KS-test may only be used to reject models, it can- 
not determine how well a model fits. Indeed, when the 
probability, Pks > 0.2, its value may not be accurate 
but the implication that the model and data are similar 
is still correct. For this reason, the relative probability 
of each model, Pr, is also evaluated directly from the 
minimisation algorithm. 



crease in the number of free parameters while not significantly 
changing the key results. 
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Table 3: Best-fit parameters for the model RLFs described in Sec. 4.3 along with their associated lu errors. N is 
the number of free parameters for each model, Smin is the minimum value of eqn. ^2|, Det(VV5') is the determinant 
of the Hessian matrix evaluated at 5 = Smin and Pr, is the relative probability of the models with respect to model 
B. A detailed description of how the errors and relative probabilities are ascribed can be found in JROO. 



4.4 Results of the parametric modelling 

Table ^ shows the best-fit parameters for the various 
models: S'min is the minimum of eqn. ^ and Pr, is the 
relative probability with respect to model B. The pa- 
rameter describing the slope of the luminosity function 
/3 is consistent with the previous results of DP90 and 
Willott et al. (2001) for the steep-spectrum RLF and 
Willott et al. (1998) for the RLF of steep-spectrum 
quasars. A 2D KS-test to the predicted L151 — z plane 
from the best-fit parameters give 0.2 < Pks < 0.4 for 
all of our models in both cosmologies suggesting good 
working models. The best-fit values for the correlation 
between ai and a2 are consistent over all models in 
the two cosmologies and are also statistically reason- 
able with a 2-D KS-test giving Pks = 0.31 (Fig. |). 
The parameters describing the Gaussian distribution 
in logj^Q D also span a narrow range between models 
and KS-probabihties of Pks ~ 0.56 (cosmology I) and 
Pks ~ 0.34 (cosmology II) suggests adequate fits to the 
data for all models. The histogram of logD / kpc and 
the fitted Gaussian arc shown Fig. ^. 

The parametric models show that model B, with a 
constant comoving space density above a critical red- 
shift is ~ 20 (~ 100) times more likely than model A, 
the model with a symmetric Gaussian decline in cos- 
mology I (cosmology II). However, one may associate 
these relative probabilities with the forced symmetry in 
model A, thus forcing a poorer fit at either low- or at 
high-redshift and these models must be considered with 



this property in mind. The best-fit parameters of model 
C are probably more relevant in discriminating between 
models with sharp cut-offs and those with constant or 
increasing space densities. 

Model C uses a 1-tailed Gaussian up to a critical 
redshift and a power-law exponent beyond this critical 
redshift, thus decoupling the low- redshift from the high- 
redshift space density. However, this model requires one 
more free parameter than models A and B. This intro- 
duces an additional factor when calculating the relative 
likelihood of this model, which is described in detail 
in JROO. 1^ The best-fit model suggests that there is 
very little evidence for an abrupt decline in the comov- 
ing space density of the most luminous low-frequcncy- 
selected sources at high redshift. Indeed the best-fit to 
the power-law exponent in cosmology I (Table ^ sug- 
gests a steady increase up to an indeterminable redshift, 
and the best-fit in cosmology II is a power-law exponent 
of '--^ thus giving a constant comoving space density, 
albeit with an uncertainty encompassing moderate de- 

§ On the assumption that most of the prior ranges cancel on 
dividing the probabihties of model A and model B, this additional 
factor is given by, 

jr= (4^)iT-^A7?, (13) 

where TV is the number of free parameters and the 11 arises from 
the number of free parameters in the comparison model (B), and 
Ari is the prior range of the additional parameter. The factor 
is of order unity and has been neglected for our analysis here. 
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Figure 4: The comoving space density <I> 
[/////pd(logioLi5i)d2;daida2d(logioi:>)] of the 
most himinous low-frequency-selected r adio sources 
using the three models described in Sec. for both 
cosmology I (top) and cosmology II (bottom). The 
dotted line is from the best-fit parameters of model 
A; the solid line for model B; and the dashed line 
represents model C. The shaded region corresponds 
to the IfT uncertainty for the power-law exponent of 
model C. This confidence region is uncertain beyond 
the limit of our dataset and thus we have excluded it at 
z > 4.41, i.e. beyond the highest redshift source in our 
samples, although as the redshift distribution in Sec. || 
shows, these uncertainties have very little influence on 
the high-redshift evolution due to the lack of available 
comoving volume. 



clines and inclines. This uncertainty on the power-law 
exponent rj is represented by the shaded region in Fig. ^ 
and also by a probability distribution function (pdf) in 
Fig. H and it is obvious that there exists a large range 
in which rj may reside. This behaviour is different from 
the pdf for the power-law exponent from the paramet- 
ric modelling of the flat-spectrum population of JROO at 
« 2 — 3cr level, where the pdf derived from the V/Vmax 
method peaked at w — 3, which is also consistent with 
their best-fit parametric models. 

The dissimilarity between the behaviour observed 
by JROO for the most luminous flat-spectrum sources 
from a high-frequency selected sample, where there was 
~ 2(7 evidence for a decline by a factor of ^ 4 between 
2.5 < z < 5, and this analysis is resolvable within the 
large uncertainties. Although the steep decline ruled 
out at the 2a level for flat-spectrum sources in JROO 
is ruled out at roughly the 4cr level here. However, if 
the difference in the best-fit power-law exponents is a 
real difference between the two radio source populations 
then it may not necessarily be due to a breakdown of 
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Figure 5: The probability density function for the 
power-law exponent rj of model C for both cosmology I 
(top) and cosmology II (bottom) , normalised such that 
the area under the curve is unity. The pdf was deter- 
mined by marginalising over the free parameters po, 
/3, Zo and zi, i.e. the four parameters with the largest 
uncertainties from the modelling (all other parameters 
were fixed at their best-fit values) . The unshaded region 
corresponds to the 90% confidence region. 

orientation-based unified schemes and we will discuss 
this in Sec. ^. 



4.5 The effect of changing the lower 
luminosity limit 

We now investigate the effect that the positioning of 
the lower limit in radio luminosity has on the outcome 
of the parametric modelling. Each model was re-fitted 
with both the lower-limit in radio luminosity increased 
by 0.25 dex and also lowered by 0.25 dex. The relative 
probabilities of the best-fit models in both cosmologies 
changed very little with this alteration. The uncertain- 
ties on the fitted parameters increased when the lumi- 
nosity limit was raised, as one might expect due to the 
small number statistics becoming increasingly impor- 
tant. The steepness of the RLF, governed by the /3 pa- 
rameter in our models was consistently « 2.1 within 
the uncertainties, suggesting that the steepness of the 
RLF in the high-luminosity regime is consistent with (3 
at lower luminosities. The power-law exponent rj for the 
high-redshift evolution in model C also varied very little 
from a best-fit value of « 0, although the large errors as- 
sociated with this parameter still make shallow declines 
and also shallow increases easily plausible. Therefore, 
we conclude that the positioning of the lower-limit in 
radio luminosity has an insignificant effect on the rela- 
tive probabilities and the fitted values of the parametric 
modelling if the modelling is restricted to only the most 
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radio luminous sources, where we are confident of com- 
pleteness in all of the samples. 

To test the susceptibility of our models to small 
number statistics we also re-fitted our models with an 
additional fake z ~ 6 source. Inclusion of the z — 6 
radio source in our models also has very little effect on 
the relative probabilities, although model A becomes 
slightly less plausible relative to model B in both cos- 
mologies as one would expect. The rj parameter is still 
consistent with a value of ?7 « and no positive evolu- 
tion in the comoving space density is required (although 
the best-fit values are rj = 0.35 and t] = 0.20 for cosmol- 
ogy I and II respectively) . The peak of the low-redshift 
Gaussian is constrained at a value of Zo ^ 2, although 
the errors on this parameter are high and correlated 
with the r] parameter. We therefore conclude that in- 
cluding a radio galaxy at z = 6 has little effect on the 
overall conclusions of our main analysis, as the small 
number statistics involved in including a single source 
at high redshift are easily compatible with the original 
best-fit models without the z = 6 source. This is also 
highlighted by the predicted redshift distribution from 
our original model C (Fig. 0), where the presence of 
a single source at z ~ 6 could easily conform to the 
predictions in both cosmologies. 



5 THE BANDED V/Vmax STATISTIC 

It is also useful to view our parametric modelling along 
with a non-parametric form in which the behaviour in 
the space density of the most luminous, high-redshift 
sources does not depend on the assumed parameterisa- 
tion. The most common method for this non-parametric 
analysis is the t^/Knax method in its banded form (Avni 

6 Bahcall 1980; see also DP90, and JROO). Unfortu- 
nately, with a filtered sample such as 6C* the actual 
significance of any value of V/14iax > 0.5 may be viewed 
as inconsequential due to the bias towards finding ob- 
jects at high redshift. However, if this banded test shows 
a significant decline in the comoving space density of 
sources with a sample biased to finding objects at high 
redshift then it would prove to be a significant detection 
of any redshift cut-off. Fig. ^ shows the banded V/Vamx 
statistic for all of the sources in the top-decade in radio 
luminosity considered in our parametric modelling. It 
is apparent that the statistic shows no significant evi- 
dence for a decline in the comoving space density out 
to z 4.5, if the highest-redshift (lensed) radio galaxy 
6C*0140+326 is included. If we now repeat the anal- 
ysis without this object there also appears to be lit- 
tle evidence of a decline, with all points lying within 
IfT of the line of no-evolution at V/Vmax = 0.5. There- 
fore, we find no evidence for a decline in the comoving 
space density of the most luminous radio sources using 
this non-parametric approach, in agreement with the 
best-fits from our parametric models. The high values 



of V/Vi-ns,x in Fig. ^ reflect a combination of genuine pos- 
itive evolution at z < 2 as the 3CRR sources dominate 
in this regime, and biases due to selection methods at 
higher redshifts where the sample becomes increasingly 
dominated by 6C*. 
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Figure 6: The banded V/Vmax statistic for the most ra- 
dio luminous steep-spectrum sources considered in this 
paper. The solid line corresponds to the statistic with 
the lensed radio galaxy 6C*0140-|-326 included in the 
analysis and the dashed line is the statistic excluding 
this source. The solid vertical lines correspond to the la 
errors on the statistic assuming a = l/\/127V (Avni & 
Bahcall 1980), where N is the number of sources above 

Zo. 



6 PREDICTED REDSHIFT 
DISTRIBUTIONS 

We now consider the redshift distribution expected 
given our best-fit values for model C. Fig. |^ shows 
the cumulative redshift distribution for our sample for 
z > 2, along with the predicted distribution for model 
C from our best-fit model. The forced symmetry in 
the estimation of these errors by fitting a symmetrical 
Gaussian to the minima in the parametric modelling 
(see JROO) means that the uncertainties may be over- 
estimated on one side and underestimated on the other. 
Therefore, although the predicted redshift distributions 
are in accord with the data the actual la errors may not 
be telling us the whole story. Thus, the pdf shown in 
Fig. 1^ is probably a better representation of the uncer- 
tainties in the parameter space we are most interested 
in. We can see that the model parameters are consis- 
tent with the data all the way through the high redshift 
regime. The upper limit on the rj parameter from our 
la uncertainty is easily consistent with the data. This 
highlights the effect of the small number statistics un- 
der consideration at the highest redshifts. The steepness 
of the RLF, combined with the sharp reduction in the 
available comoving volume due to the flux-density limits 
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of the samples means that, even without any dechne in 
the comoving space density of radio sources, the number 
of sources one would expect to find is < 1. 




Redshift 

Figure 7: Cumulative redshift distributions at z > 2 
for model C in cosmology I (top) and cosmology II 
(bottom). The histogram represents the actual cumu- 
lative redshift distribution of the sources in the 3 sam- 
ples used, the hashed regions are the VN errors on the 
binned data. The solid curve is the predicted distribu- 
tion from the best-fit values of model C. The two dot- 
ted curves represent the curves for the la uncertainty 
on the liigh-rcdshift exponent r], using the uncertainties 
generated in our parametric modelling. 



7 THE RLF OF STEEP- AND 

FLAT-SPECTRUM RADIO SOURCES IN 
THE CONTEXT OF UNIFIED SCHEMES 

Our work on the radio luminosity function of both 
steep-spectrum radio sources (this paper) and flat- 
spectrum radio sources (JROO) has highlighted the dif- 
ficulty in establishing the comoving space density of 
these sources at redshifts which push to the limits of 
flux-density-limited radio samples. Although the uncer- 
tainties on the fitted parameters describing the high- 
rcdshift behaviour of both these types of radio source 
arc quite large they are markedly different in their mean 
value for the high- redshift power-law exponent. Indeed 
the evidence from our parametric modelling favours a 
steady decline for the flat-spectrum population by a fac- 
tor ~ 4 between 2.5 < z < 5 (JROO), whereas analysis 



of the low-frequency samples in this paper lead us to a 
best-fit model with a roughly constant comoving space 
beyond a peak redshift around z w 2.5. 

This discrepancy is just resolvable within the un- 
certainties of both models, so it is plausibly a statistical 
fluke that we arc seeing different evolutionary behaviour 
between the two populations of radio sources. However, 
other plausible explanations exist and should be inves- 
tigated further. 

The most obvious difference between the flat- and 
steep-spectrum samples is the selection frequency. The 
high selection frequency (2.7 GHz) for the flat-spectrum 
sample means that the majority of objects are core- 
dominated Doppler-boosted quasars and, as argued by 
JROO, Giga-Hcrtz Peaked spectrum (GPS) sources may 
also make an important contribution. If a signiflcant 
fraction of the most luminous flat-spectrum sources are 
indeed GPS sources, then at an observed frequency of 
2.7 GHz, the radio spectrum of a GPS source may have 
steepened by an amount sufficient for it to drop below 
the sample flux-density limit or for it to have been re- 
jected from the sample on account of its steep observed 
spectral index {a > 0.5) if it is at high redshift. Indeed 
even if the core-dominated, Doppler-boosted sources arc 
a dominant part of the top-decade of the flat-spectrum 
luminosity function, and assuming a conservative esti- 
mate for the fraction of GPS sources, a decline in the 
comoving space density would be observed at high red- 
shift. This is because of the drop off in GPS sources 
due to K-correction effects, allowing a roughly constant 
comoving population of Doppler-boosted sources. This 
point was made by JROO and they suggest a model for 
the RLF is needed in which both the Doppler-boosted 
and GPS populations are explicitly considered. The 
banded V/Vmax method used by JROO was able to ac- 
count for the spectral shape of individual sources and 
although small number statistics make it very difficult 
to be certain (see the broad 90% confidence region in 
their Fig. 9b), the presence of some high- redshift decline 
in one or more of the Doppler-boosted and GPS pop- 
ulations appears robust to spectral-dependent selection 
effects. 

Another difference between the flat- and steep- 
spectrum sources is the steepness of the luminosity func- 
tion in the top decade. JROO found values in the range 
1.2 < (3 < 1.75 for the top decade of the flat-spectrum 
population, whereas values of /3 > 2.0 have been de- 
termined for the top decade of steep-spectrum RLF in 
this analysis and over a broader range in luminosities 
in previous work (e.g. Willott et al. 1998; Willott et al. 
2001). DP90 also found a marginally steeper luminosity 
function for their steep-spectrum population (/3 ~ 2.2) 
compared to their flat-spectrum sample (/? ^ 2.0). 

The difference in both the high-redshift decline 
and values of /? between the flat- and steep-spectrum 
sources can be resolved by considering the concave 
shape of the RLF (Wfilott et al. 2001). The most lumi- 
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nous flat-spectrum population probes a different part 
of the underlying RLF than the most luminous low- 
frequency-selected sources (JROO). Thus, if the posi- 
tion of high-frequency selected flat-spectrum sources on 
the underlying RLF are at luminosities which are lower 
than those of the most luminous low-frequency-selectcd 
sources and lie on or near the 'break' in the RLF at 
logio-^151 ^ 26.5 for cosmology I (Willott et al. 2001), 
then this automatically supports the steeper values de- 
rived for P from our steep-spectrum samples. In other 
words, the slope in the RLF steepens with increasing 
intrinsic luminosity, so either Doppler boosted or GPS 
sources in high-frequency samples will be drawn from 
a population where the slope in the underlying RLF is 
flatter. Indeed, it seems likely that the Doppler boosted 
and GPS objects are drawn from different points on 
the underlying RLF, so there may well be differential ef- 
fects: JROO for example, have suggested this as the cause 
of an increasing fraction of GPS sources with increas- 
ing luminosity in high-frequency selected, flat-spectrum 
samples. 

The difference in the evolution of high-redshift 
comoving space density between the flat- and steep- 
spectrum populations can now be imderstood if there is 
a decline in the comoving space density at high redshift 
for lower-luminosity low-frequency sources. This would 
be very important to confirm because the luminosity 
density (/ Lpdlogj^Q-L) of the population is dominated 
by the space density near the break of the RLF. It is 
impossible to probe this break region directly with the 
complete flux-density-limited samples currently avail- 
able. 

We therefore conclude that, although the high- 
redshift space density and the steepness of the RLF may 
be different for flat- and steep-spectrum sources in the 
top-decade of luminosity, this can easily be explained 
by some combination of two effects: flrst, a mixed pop- 
ulation of Doppler boosted and GPS sources in the 
high-frequency selected flat-spectrum samples; second, 
and arguably more importantly, the flat- and steep- 
spectrum populations lie on different parts of an under- 
lying RLF which is concave (Willott et al. 2001). Thus, 
orientation-based unified schemes are not challenged by 
the comparison of constraints on the high-redshift evo- 
lution of flat- and steep-spectrum radio sources. 



8 CONCLUSIONS 

The main conclusions we have drawn from our para- 
metric modelling analysis of the most luminous low- 
frequency-selected radio sources are: 

• Their low-redshift evolution follows the same track 
as that of fiat-spectrum radio sources, optically selected 
AGN (e.g. Schmidt, Schneider & Gunn 1995) and X-ray 
selected AGN (e.g. Miyaji, Hasinger & Schmidt 2000), 



with an increase in comoving space density by a factor 
>100 over 0<z< 2.5. 

• Their steep increase in space density slows down at 
z > 2.5 in spatially-flat cosmologies with either f^A = 
or Qa = 0.7. 

• The saturated Gaussian model B with a constant 
comoving space density beyond a peak redshift, is « 25 
(w 100) times more likely than model A with a Gaus- 
sian rise and cut-off for cosmology I (II) . Although this 
ratio may partly reflect the forcing of the low-redshift 
increase to be symmetric with the high-redshift decline 
in the comoving space density of the most luminous ra- 
dio sources, it does illustrate the signiflcant difference 
between radio-loud AGN and optically selected quasars 
[which do have a roughly Gaussian evolutionary be- 
haviour, e.g. Schmidt et al. (1995); Kennefick, Djorgov- 
ski & de Carvalho (1995)]. 

• The best-fit model is consistent with a constant 
comoving space density beyond the peak ai z 2.5 to 
an indeterminable redshift, although the uncertainties 
span a wide range covering both moderate declines and 
continuing shallow inclines. 

• The difference between the high-redshift space den- 
sity for fiat- and steep-spectrum sources is marginally 
significant, although there remain worries concerning 
the mixture of Doppler Boosted and GPS sources in 
the flat-spectrum population. An explanation exists 
in which the high-frequency flat-spectrum sources are 
probing the underlying RLF at fainter intrinsic lumi- 
nosities than the top-decade discussed in this paper. 
This could also resolve the difference in (3 between 
the two populations, because at lower-luminosities the 
sources are edging towards the break in the luminosity 
function where j3 becomes flatter (Willott et al. 2001). 
This is important to conflrm since the tentative evi- 
dence for a decline in the flat-spectrum population sup- 
plies indirect evidence of a decline in the luminosity 
density in the radio source population at high redshift. 
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